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Abstract 

This paper presents a general method for incorporating state 
transformation in event driven explicit simulation. One inher- 
ent assumption in this type of simulation algorithm is the state 
independence which allows the algorithm to process the 
states independently in an event driven manner at the transis- 
tor level Numerical problems arise when an inappropriate 
state representation of the circuit, in which the states are not 
truly independent, b chosen. In principle, any similarity 
transformation of the state equation can be employed to 
transform the circuit into a more convenient state space for 
numerical solution. This paper develops a systematic scheme 
to derive an appropriate state transformation, and to incorpo- 
rate the state transformation in such a way to maintain the 
efficiency of event driven explicit simulation algorithms. 

1. Introduction 

With design trends continuously moving towards higher clock 
speed, higher density, and higher levels of integration on chip, the 
need for fast circuit simulation of large and complex integrated cir- 
cuits is well known. There have been two major approaches in 
improving the efficiency of (circuit) transient simulation. Common 
to all approaches is the use of circuit partitioning and event driven 
simulation to exploit the spatial sparsity and the temporal latency of 
large circuits. One approach is to simplify circuit topology to a set of 
primitives, especially for MOS circuits, in addition to employing 
simplified device models, such as piecewise quadratic models, to 
derive the analytical solutions of the node voltages. This approach is 
embodied in the simulator ILLIADS [3]. The other major approach 
employs simplified device models to speed up the computation of the 
numerical solutions during a transient analysis. Piecewise constant 
(PWC) device models are employed in SPECS [4] and ADAPTS [5]. 
More general PWL models have been used in various timing analyz- 
ers and simulators [1][6]. The Adaptively Controlled Explicit Simu- 
lation (ACES) algorithm [2], has been developed to exploit the 
efficiency of explicit integration algorithm white ensuring the 
numerical stability of the integration at an incremental cost instead 
of the full cost of an implicit algorithm. The basic idea is to incre- 
mentally approach the steady state while ensuring the stability and 
efficiency of the integration algorithm. 

In event driven simulation, an implicit assumption is that the 
states are independent so that they can be processed separately. If 
some of the states are strongly coupled or dependent, these states 
have to be processed as a group. In other words, the choice of the 
state variables for a state representation of a given circuit is not 
appropriate in the sense that the chosen states are not truly indepen- 
dent. This situation arises quite often in extracted circuits, which 
contain numerous parasitic resistors of small values. For example. 



given a resistor as shown in Fig. 1 , the two capacitors at the two 
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Fig. 1; A simple illustration of the suite dependency problem asso- 
ciated with a small value resistor 



nodes of the resistor are usually chosen to be the state variables in a 
typical circuit formulation such as Modified Nodal Analysts (MNA). 
In the degenerate case of a zero valued resistor or voltage source, the 
two states are actually the same. Even if the value of the resistance 
is very small, i.e. the two nodes close to being shorted, then the two 
chosen state variables are not truly independent In the following, a 
number of numerical problems associated with the case of small 
resistors mentioned above are discussed. 

Matrix ill-conditioning 

The first one is the numerical ill conditioning problem due to the 
resistor being stamped into the MNA circuit matrix equation as a 
conductance term with G = 1 /R in the admittance form 

i R = Gv R = C(v I -v 2 ). 

The reason that resistors are stamped as conductors is to avoid 
introducing extra unknowns into the MNA matrix equation, in which 
the capacitor voltages and inductor currents are chosen as state vari- 
ables. However, if R is small, then G is large and it may degrade the 
numerical conditioning of the matrix equation. One way to avoid this 
problem is to stamp the resistor in its impedance form as 



v R -Ri R = v r 



-Rir 



Note that this impedance form introduces the extra unknown cur- 
rent variable i R into the matrix equation. Note further that this 
stamping does not address the state dependency problem though it 
may alleviate the numerical conditioning problem. 

State Dependency in Explicit Simulation 

The state dependency problem is evident in explicit integration 
algorithm such as Forward Euler (FE). Given the Branch Constitu- 
tive Relation (BCR) of a capacitor as i c = Cv c . the FE stamp of 
the capacitor at each time / is given by 



where r - Ar is the previous time point and Ar is the time step. The 
above equation represents the BCR of a voltage source, and the 
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equivalent circuit of the simple structure shown in Fig. I using FE 
integration is shown in Fig. 2 below. As illustrated by Fig. 2, if R is 




v C2 {t-At) 



Fig. 2;The equivalent circuit of the circuit shown in Fig. I at time t 
using an explicit (Forward Euler) integration for the capacitors. 



ulation (ACES) algorithm as an implementation vehicle. This paper 
is organized as follows. In section 2 of the paper, a simple example 
is presented to review the basic ACES algorithm as an implementa- 
tion vehicle and illustrate the numerical problems due to state depen- 
dency. The main results of the paper are presented in section 3. This 
section discusses a general and systematic method to incorporate 
state transformation into event driven explicit simulation, especially 
its implementation in ACES. Section 4 presents some additional 
examples and simulation results on real circuits using state transfor- 
mation. Section 5 summarizes the paper. 



small or zero in the degenerate case, the equivalent circuit includes an 
illegal loop of voltage sources. This state dependency problem will 
affect any explicit integration algorithm that simply update the capac- 
itor voltages at a given time using only values at previous time points. 

State Dependency in Implicit Simulation 

The problem of state dependency is handled automatically in an 
implicit integration scheme because all the states are 'Implicitly" pro- 
cessed together at the same time. For example, using Backward Euler 
(BE) as the implicit integration scheme, the stencil of a capacitor in 
the MNA matrix equation at each time / can be written as 



v c (0 = ^ c ('> + w c (f-A/) 



In circuit theoretic terms, this equation represents an equivalent 
circuit consisting of a resistor with conductance G eq - (A/)/C in 
series with a voltage source V = v c (t -At) . Using these equiva- 
lent circuits for the two capacitors, the circuit shown in Fig. 1 can be 
redrawn as shown in Fig. 3 below. In this case, even if R = 0 , the 
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Fig. 3:The equivalent circuit of the circuit shown in Fig. 1 at time t 
using an implicit (Backward Euler) integration for the capacitors 

two equivalent voltage sources are still decoupled by the two equiva- 
lent conductors, and these voltage sources do not form an illegal loop. 

The fundamental problem is then how to handle this state depen- 
dency problem without resorting to implicit integration algorithms 
for efficiency reasons. One simple solution is to preprocess a given 
circuit and delete from the circuit description all the resistors with 
small values below a specified tolerance. However, this approach may 
require setting different tolerances fordifferent circuit configurations. 
Moreover, in some cases, it may be necessary to retain these small 
resistors in the circuit description. For example, given a long RC 
interconnect with many taps along the interconnect, it is necessary to 
retain the small resistors in each RC section of the interconnect to get 
the total resistance loading of the line. This problem obviates the need 
to develop an algorithmic approach, other than circuit pruning, to 
handle this problem. It is the goal of this paper to present such an 
approach using the event driven Adaptively Controlled Explicit Sim- 



2. State Dependency in Event Driven Explicit Simulation 

Consider a simple 2RC circuit as shown in Fig. 4 This simple exam- 

R1 CD 5\ O 
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Fig. 4; A simple 2RC example to illustrate the state dependency 
problem 

pie is used to review the basic ACES algorithm and to illustrate the 
problem of state dependence. The state equation of the circuit can be 
written as 



C, 0 
0 C, 



~G?. + *J R 



(i) 



where the capacitor voltages have been chosen to be the state vari- 
ables. Let V in be a unit step function, C, = C 2 = IF,/?, = lfi, 
and R 2 = 10 Q . Then the state equation can be written as 



'ii = T-o+io 6 ) io 6 i| v ii. 



V.- 



(2) 



Assume that the circuit is initially at rest, i.e v, = v 2 = 0 at 
t = 0 . Also assume a unit step input, i. e. V in - 1 and V in = 0 
at t = 0 . Then, applying the ACES algorithm to compute the time 
to quiescence [2] for the two states yields 
Ar j s -(v 1 (0))/(v,(0))« 10 -6 and A/ 2 = 0 , which implies that 
the input excitation has not reached the second node at this time 
point. In other words, one is interested onlv in non-zero time step at 
each time point. Hence, A/ = A/j = 10 

Continue the ACES algorithm to the next time point, 
t x - t 0 + At = 10" 6 . Note that v, is now in quiescence and it 
should not dictate the computation of the next time step. Hence the 
time to quiescence for v 2 is given by 

A/ = -(v^,))/^/,)) = 1 + 10-*. 

Then, at the next time point, / 2 = t x + At , both states will be qui- 
escence and the simulation is completed. The result of this simple 
simulation for v 2 is shown in Fig. 5 together with the result obtained 
by AS/X, the IBM internal circuit simulator. Fig. 5 shows the result 
obtained by ACES with a voltage error tolerance e v = 0.001 , where 
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Fig. 5: ACES simulation result (without state transformation) as 
compared to AS/X for the simple 2RC circuit shown in Fig. 4 with 
an errpr tolerance ofO.OOL 

e v = -y(At) 2 . Even with this fairly tight tolerance, the ACES result 
shows "significant discrepancy as compared to the AS/X result. This 
problem is due to the fact that the two states are treated as indepen- 
dent and processed one at a time. 

From the circuit point of view, when enters quiescence, ACES 
replaces C, by a zero valued current source, i. e. C, becomes an 
open circuit and the time constant of this new circuit is given by 
t = +R 2 )C 2 = 1 + 10" 6 . This is exactly the time to quies- 
cence or the time to steady state for this new circuit, assuming a first 
order approximation for v 2 . However, this is not the correct time con- 
stant for the original circuit. Ignoring the effect of the small resistor 
/?2 . the correct time constant is given by t = /?|(Cj + C 2 ) = 2. 
This discrepancy illustrates the problem in assuming the indepen- 
dence of strongly coupled state variables. 

In order to resolve this problem, an appropriate state space repre- 
sentation of the circuit must be chosen in such a way that the states of 
this space are truly independent. One obvious choice is the eigen- 
space of the circuit. Compute the eigcn -decomposition of the system 
matrix as 



A = 



1 1 



-U + 10 6 ) 10 6 
10 6 -10 6 



>tat- 



(3) 



where 



T = 



1 1 




i i 


75 Jl 


.r' = 




1 1 




1 1 


.72 72 




.72 72. 



,and A = 



-2(10 6 ) 0 
0 -0.5 



(4) 



Define the transformed state variables as 
\y \ ^2] = [ V 1 v 2] ^ * ^ e trans f° nne d state equation can then 



be written as 





= A 
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Now apply the ACES algorithm to this system. At / = 0, 
y x s y 2 = 0 . From the state equations . it can be seen that y ( * 0 
j?2 * 0 , y, * 0 , and y 2 * 0 . Hence the time step can be computed as 



(5) 



Af = min(At At*,) = min\ 



At the next time point /, = A/ = O.SCIO" 6 ), y x = y 2 = 0. 
Hence y, will stay in quiescence and the time to quiescence for y 2 
can be computed as 

A/ = ~ = 0.5 

In this case, the ACES algorithm captures the correct time con- 
stant of the original circuit. This simple example has clearly illus- 
trated the need to select an appropriate state space in which the state 
variables are truly independent so that an event driven explicit algo- 
rithm such as ACES can be applied to provide accurate simulation 
results. 

3. State Transformation for Event Driven Explicit Simula- 
tion 

The ei gen- space representation is an ideal space for event driven 
explicit simulation because the modes of the circuit are decoupled 
and can be processed sequentially. However, this approach is clearly 
impractical due to a number of reasons. First the computation of the 
eigen-decomposition is expensive. Moreover, once the eigen- 
decomposition is known, the solution can be computed directly with- 
out resorting to numerical integration. Avoiding the eigen-decompo- 
sition is exactly the goal of numerical integration algorithm. So the 
question is how to determine an appropriate state transformation. 
One possible approach is to perform an approximate or inexact 
eigen-computation. This alleviates the cost of a full eigen-computa- 
tion. However, for non-linear circuits, an eigen-computation needs to 
be performed for every new linearization of the nonlinear equations. 
Another approach is to perform a partial eigen-computation. This is 
a basically a model reduction problem, especially for large linear cir- 
cuits. This approach is a standard technique for reduced order mod- 
eling of large linear circuits in general nonlinear simulation of 
circuits with large linear subcircuits. The advantage of this approach 
is that it can prune out the non dominant eigen-states, such as the 
state due to the large eigenvalue (2(10 6 )) in the 2RC example 
above. However, this approach is not practical for nonlinear circuits. 

In general, there is a need to have a simple transformation matrix 
in order not to degrade the efficiency of event driven explicit simula- 
tion technique such as ACES. Note again that one key assumption of 
circuit level (not logic level or macromodel or block level) event 
driven simulation is that the states are independent or DC decoupled 
so that the excitation can be propagated from the input source to the 
various nodes in the circuit in an event driven manner. For example, 
the unit step excitation from the input source of the above 2RC cir- 
cuit will propagate first to node 1 , and then to node 2. However, when 
the two nodes are strongly DC coupled, the algorithm will yield erro- 
neous results as discussed in the previous section. The key idea is to 
process the two states simultaneously. Note that the eigen-transfor- 
mation described in the previous section transforms the input matrix 
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[l o] T 10 [l/-/2 \/j2^' m tne ori S inal slale space* lne input 

directly affects only v } . In the transformed state space, the input 

directly affects both y, and y 2 simultaneously. Hence the key idea 

in constructing a simple transformation matrix is to make sure that all 
the transformed states will be dependent on all the original states. In 
this way, any input that affects at least one state in the original space 
will affect all the states in the transformed space. In this way, the set 
of coupled transformed states can be processed simultaneously dur- 
ing the "electrical" event propagation. For the 2RC circuit example, 
one such simple transformation matrix is given as 



•[:-:] 



(6) 



This matrix looks like a simple scaling of the transformation 

matrix 7" 1 (the inverse eigenvector matrix) discussed in the previous 
section. In general, this is not the case. The transformation matrix in 
Eq. (6) is simply to make every transformed states dependent upon all 
the original state variables. In order to bring out the effect of the 
capacitor values, let *s choose C, =0.5 and C 2 = 2 . Also select 

the capacitor charges, denoted by q x and g 2 as the state variables in 

the original space, i. e. [4, q^ T = [CjV, C 2 v^] r . And the trans- 
formed state variables are defined as 



■LrW = rHh 

2 J [fcJ L« UL?2 



An ACES simulation is then performed for this circuit. Note that 
the times to quiescence (an "electrical" event in ACES) are computed 
for the transformed state variables. The local truncation error control 
is performed for the original state variables (the node voltages) in this 
case. The waveform at node 2 as obtained by the ACES simulation 
with state transformation is shown in Fig. 6 together with the AS/X 

VI (V) 




Fig. 6: ACES simulation result ( with state transformaxion) as com- 
pared to AS/X for the simple 2RC circuit shown in Fig. 4 

simulation result Note that the state transformation helps ensure the 
accuracy of ACES simulation. 



In general, if there are no coupling among the states, then the 
transformation matrix is simply an identity matrix. Suppose there are 
n coupled states, then the submatrix corresponding to these n states 
will be of the form 



+1 -J -J ... _i 
+1 +1 -1 ... -1 

T nxn = +1 +1 +1 ... -I 



+1 +1 +1 ... +1 



For example, given a circuit with 7 states among which there are 
two sets of coupled states: {2,3,4} and (6,7), then the transformation 
matrix for this circuit can be written as 



10 0 00 0 
0 1-1-10 0 
0 1 1-10 0 
1 10 0 
0 0 10 
0 0 0 1 -1 
0 00 1 1 



0 1 
0 0 
0 0 
0 0 



This simple structure allows the transformation matrix to be con- 
structed in a systematic manner during the processing of the circuit. 
For example, a simple two phase procedure to construct this type of 
transformation matrix can be described as follows. First initialize the 
transformation matrix T to be an identity matrix. Then in the first 
phase of processing the circuit, check if the terminals (which are 
assumed to be states for simplicity) or nodes of a given resistive ele- 
ment (or device) are coupled. If the nodes are coupled, then enter the 
appropriate nonzero entries, which corresponds to the indices of the 
nodes in the state vector, in the transformation matrix. The purpose 
of this phase is to record the coupling information in the matrix. In 
the second phase, the transformation matrix will be processed row by 
row to determine different groups of coupled states. Then the sub- 
matrices of the type described above for these groups of coupled 
states will be recorded in the transformation matrix. 

The next question is how to incorporate state transformation into 
ACES. A straightforward approach is simply transform the original 
state equation into a new state equation and perform the analysis in 
the transformed space. However, it should be noted that the trans- 
formed space is used to compute the time to quiescence for the trans- 
formed states while the local truncation error is still performed in the 
original space. Therefore, it is simpler to impose the constitution 
relation or the quiescence condition for the transformed states in 
terms of the original states to solve for the necessary variables in the 
original space. Then the variables in the transformed space can be 
computed by applying the transformation matrix to the original 
states. These ideas will be illustrated below using a simple Modified 
Nodal Analysis (MNA) circuit formulation in the derivative space 
[2]. 

In order to simplify the presentation and to bring out the differ- 
ence between the basic ACES algorithm and the modified algorithm 
to handle state transformation, assume that there are no inductors. 
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current sources, floating capacitors, and floating voltage sources. 
These circuit elements can be incorporated into the formulation in a 
straightforward manner once the basic algorithm has been estab- 
lished. Under these assumptions, during a quiescence period [2], the 
MNA formulation of the circuit in the derivative space can be written 
as 



Gcc O cv E cc 0 
G VC G VV 0 E y 

0 U u 0 0 









0 








0 












h 







(7) 



where G cc , G i 



G vc , and G vv are submatrices of the conductance 
matrix formed by resistive elements; E cc , E qc , E v , and U v are sub- 
matrices of Vs and O's to pick out the appropriate capacitor/source 
voltages or currents; K c and h c represent the changeable portion of 
the MNA equation during the simulation depending upon whether the 
states are in quiescence or not. In general, the first two rows of Eq. (7) 
represents Kirchoff Current Law (KCL) equations at all voltage and 
capacitance nodes. The third row represents the Branch Constitutive 
Relations (BCR) and/or the quiescence conditions for the capacitors. 
The fourth row represents the BCRs for the voltage sources. 

Now define the state transformation y c = Tq e . Initially when all 
the states are non-quiescent, the BCRs for the transformed states can 
be written as y c = TCv c . In circuit theoretic terms, the transformed 
states can be considered as voltage controlled current sources. In this 
case, the stamps in the MNA matrix of Eq.(7) can be set as follows: 
K c = TC , E qc = 0 , and h c = j? c . Suppose at time / , the j th 
(transformed) state enters quiescence. Then y c j(t) ~ fj# c (0 = 0 
and the stamps in the MNA matrix can be modified as kL = 0, 
e* c j * t j , and h cJ = 0 , where tj is the j th row of T , kf- is the 
j th row of K c , and e* c ■ is the j th row of E qc . In this case, the 
transformed state in quiescence can be treated as a zero valued cur- 
rent controlled current source. In this case, partition the vector of 
transformed states as 

















Jk 



(8) 



where y cn denotes the subvector of transformed states that are not in 
quiescence, and y c$ denotes the subvector of transformed states that 
are in quiescence; T n and T s are appropriate partitions of the trans- 
formation matrix T. Then the MNA equation as shown in £q.(7) can 
be rewritten as 



Gee G c 



0 
0 



0 
0 

T, 
0 



(9) 



where the third row of Eq.(9) describes the non-quiescent trans- 
formed states, and the fourth row describes the quiescent transformed 
states. The modified ACES algorithm with state transformation is 
summarized in the pseudocode below. 



Algorithm with state transformation : 

1. Update sources 

2. Update original state variables (for all states): 

q e (t-) s ? c (/-A/) + ? c (/-Ai)A/ and 
v c (f-) = v c (f-A/) + v c (f-Af)Af 

2b. Update the transformed states: 

)> c (/-) = Tq c {t~) for quiescent states, and y c (/-) = 0 for 
quiescent states 

3. Check quiescence conditions for y c 

4. Setup the MNA matrix and solve for v c (/+) , q c {t) , and 
WO- 

5. Compute q c (t+) = Cv c (/+) andv c (t) = C' X q c {t) for all 
states 

5b Compute y c (t+) = Tq c {t+) and y c {t) = Tq c {t) for 
non-quiescent states 

6. Compute time to quiescence for non-quiescent states, 

= -W>/<y c (/))i 

7. Set A/ = min{t ) 

In the above algorithm, the quiescence conditions are checked for 
the transformed states, not the original states. The advantage of this 
algorithm is that there is no need to explicitly construct the trans- 
formed space. The solution procedure is still performed in the origi- 
nal space with minimal modification of the matrix equation to 
impose the BCRs and/or the quiescence conditions for the trans- 
formed states. In the next section, some examples will be presented 
to assess the overhead of the modified ACES algorithm and its accu- 
racy as compared with the basic ACES algorithm. 

4. Examples 

The first example is a simple nonlinear circuit consisting of a 
driver driving four buffers. Each connection between the driver and 
a given buffer is an 8 RC interconnect model. The schematic of the 
circuit is shown in Fig. 7. This example is to illustrate the complies- 

zm 
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Fig, 7; A simple 5 inverter circuit to illustrate the complication due 
to circuit partitioning when using state transformation 

tion due to circuit partitioning which may couple a port to an internal 
state of a given partition. This situation arises when parasitic resis- 
tances are added to the terminals of a MOSFET model. For example, 
there are five partitions PI to P5 in the circuit shown in Fig. 7. Ports 
numbered 1 to 3 will be coupled to the internal states of partition PI 
if the resistor at the end of the lines Tl to T3 are shorted. In a stan- 
dard ACES algorithm, a partition is processed only if a port attached 
to the partition has an event A port event is initialed when the volt- 
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age slope at the port has an appreciable change. Otherwise, the parti- 
tion is not processed even if the port information is updated when 
another partition connected to the same port has an event. For exam- 
ple, assume partition P2 has an event at a given time point. In process- 
ing partition P2, the voltage and the slope at port 1 will be updated. If 
the slope has a significant change (above a given tolerance), then par- 
tition PI will be processed at the same event time. Otherwise, PI will 
not be processed. However, in the case port 1 is coupled to an internal 
state of PI , PI will have to be processed regardless of whether PI has 
a significant change in voltage slope or not. This processing is neces- 
sary to take the coupling into account when computing all the neces- 
sary information for the next event. A sample run is performed for 
this circuit with the first, last, and two middle resistors of the lines are 
shorted to produce different modes of coupling among the internal 
states and the ports of PI. In this case, the ACES simulation without 
state transformation requires 30971 8 events and 34.36 CPU seconds. 
The ACES simulation with state transformation requires only 3822 
events and 2.77 CPU seconds. The ACES simulation results are com- 
pared with AS/X simulations and shown in Fig. 8 at the output of par- 
tition P2. 
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Fig. 8; ACES simulation results (with and without state transfor- 
mation) as compared to AS/X at the output of partition P2 as 
shown in Fig. 7. 

The last example is an industrial design of a 4 bit ALU. The circuit 
consists of 282 MOSFETs and 156 nodes. Without parasitic resis- 
tances added to the MOSFET models, the ACES simulation requires 
only 7298 events and 1.17 CPU seconds while the AS/X simulation 
requires 50 CPU seconds. With a number of parasitic resistances 
added to the circuit, the AS/X simulation requires about the same 50 
CPU seconds. In this case, the ACES simulation with state transfor- 
mation requires 1 7640 events and 5.73 CPU seconds. The waveforms 
at an output node of the circuit as obtained by ACES simulation with 
and without parasitic resistances and by AS/X are shown in Fig. 9. 
Note that the difference between the ACES waveforms is due to the 
state transformation. In the case of no parasitic resistances (and no 
transformation), the time to quiescence is computed for the original 
state. In the case of state transformation, the time to quiescence is 
computed for the transformed states. 

5*Summary 

This paper has identified the dependency of the states in the state 
space description of a given circuit as the cause for some numerical 
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Fig. 9: ACES simulation results (with and without parasitic 
resistances) as compared to AS/X at the output of an industrial 
4 bit ALU circuit, 
problems in event driven explicit simulation algorithms such as 
ACES. State transformation has also been proposed to address this 
problem. A novel algorithm has been developed to derive systemat- 
ically a simple state transformation and to include such a transforma- 
tion in an efficient manner. A number of examples have been 
included to illustrate the concept and the complication in using state 
transformation as well as to demonstrate its effectiveness. 
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